Classical and Quantum Chaotic Behaviors of Two Colliding 

Harmonic Oscillators 



Qing-Rong Zheng 1,2 , Gang Su 1,3 and De-Hai Zhang 1,4 
Physics Department, Graduate School, Chinese Academy of Sciences, P. 0. Box 3908, Beijing 

100039, China 

2 Department of Physics, The University of Hong Kong, Pokfulam Road, Hong Kong 
3 Institut fur Theoretische Physik, Universitdt zu Koln, Zulpicher Strasse 77, D-50937 Koln, 

Germany 

Center of Theoretical Physics, CCAST (World Laboratory), P.O. Box 8730, Beijing 100080, 

China 

Abstract 

We have systematically studied both classical and quantum chaotic behaviors 
of two colliding harmonic oscillators. The classical case falls in Kolmogorov- 
Arnold-Moser class. It is shown that there exists an energy threshold, above 
which the system becomes nonintegrable. For some values of the initial en- 
ergy near the threshold, we have found that the ratio of frequencies of the two 
oscillators affects the Poincare sections significantly The largest Lyapunov 
character exponent depends linearly on the ratio of frequencies of the two os- 
cillators away from the energy threshold in some chaotic regions, which shows 
that the chaotic behaviors of the system are mainly determined by the ratio. 
In the quantum case, for certain parameters, the distribution of the energy 
level spacings also varies with the ratio of frequencies of the two oscillators. 
The relation between the energy spectra and the ratio of frequencies of the 
two oscillators, the interaction constant, and the semi-classical quantization 
constant, is also investigated respectively. 



I. INTRODUCTION 



Many simple systems show rather complicated classical and quantum chaotic behaviors 
due to collision. A well-known example is the billiard system which describes the collision 
of a free mass point with various shapes of boundaries[l]. Another example is the occur- 
rence of the chaotic scattering due to collisions[2]. We know that there are many processes 
related to the collision phenomena in the real world, such as chemical reactions [3], celestial 
mechanics[4], charged particles in an electromagnetic field[5], processes in hydrodynamics [6], 
scattering in atomic physics[7], and transportation problems in mesoscopic physics [8], etc. 
A lot of papers are devoted to investigate collision problems in literature. The study of 
colliding phenomena, however, is far from complete, even for the elastic collision. 

The elastic collision can be assumed to be described by a so-called point contact interac- 
tion or ^-interaction. In the regime of quantum mechanics, several one-dimensional systems 
with 5 -interaction potential are exactly solved[9, 10, 11]. Especially, the solutions of many 
body systems with (^-interaction have deep influences on physics even on mathematics[10, 11]. 
A lot of people have studied the following Hamiltonian 

ff = E^ + EV(^-4 (!) 

i ij 

with Pi the momentum of the ith particle. A simple generalization of eq.(l) is that each 
particle has its own potential V(x) such that 



This Hamiltonian system, in general, can be exactly solved neither in classical mechanics 
nor in quantum mechanics. The system of two colliding harmonic oscillators is just the 
nontrivial reduction of the Hamiltonian eq.(2). 

In this paper, we will study systematically, both in the sense of classical mechanics and 
quantum mechanics, the simplest nontrivial reduction of the extended Hamiltonian eq.(2): 
two colliding harmonic oscillators, i.e., two harmonic oscillators with ^-interaction. This sys- 
tem is conserved, and the classical case falls in Kolmogorov-Arnold-Moser(KAM) class[12]. 
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Besides, there exists an energy threshold, above which the system becomes nonintegrable. 
For some values of the initial energy, we have found that the ratio of frequencies of the 
two oscillators affects the Poincare sections significantly. The largest Lyapunov character 
exponent (LCE) depends linearly on the ratio of frequencies of the two oscillators away from 
the energy threshold in some chaotic regions. In the quantum case, we have studied the 
distribution of the energy level spacings and the density of probability determined by the 
wave function, and the relation between the energy spectra and the ratio of frequencies of 
the two oscillators, the interaction constant and the semi-classical quantization constant, 
respectively. The paper is organized as follows: The classical case of the two colliding os- 
cillators will be discussed in Sec. II, and the quantum case will be studied in detail in Sec. 
III. The discussions and concluding remarks will be given finally. 



II. THE CLASSICAL CASE 

As stated above, we will study the two colliding oscillators, the simplest nontrivial re- 
duction of the Hamiltonian eq.(2). The Hamiltonian of the system can be written as 

H = H 1 + H 2 + H', (3) 

where 

H ^^ + ¥^ +b)1 - (4) 

and H' is the point contact interaction 

H' = VS{x 1 -x 2 ), (6) 

where k±, k 2 > are stiffness constants, 2b is the distance between two oscillators, and V is 
the interacting constant. 

The classical dynamic equations of the system can be in principle written down, but the 
resulting equations contain the first order derivative of 5-function explicitly. Since numerical 



methods can not deal with such an indefinite quantity directly, an approximate expression 
of 5-function must be used in the situation. Whereas in the quantum case it is not so, i.e., 
the 5-function can be treated exactly. Based on this consideration and for the consistency of 
the two cases, we will use the collision condition to replace the (^-interaction in classical case. 
Such a replacement is equivalent to the case that the (^-function is directly used. Therefore, 
the equations of motion of our system are written down as follows 

xi = Pi/m 1 

A = -k^ + b) ^ 

x 2 = P 2 /m 2 

P 2 = -k 2 (x 2 - b). 

for Xi 7^ x 2 , and for X\ = x 2 with an additional condition: the exchange of Pi and P 2 at X\ = 
x 2 . In this way, it is easy to understand that the behaviors of the system are independent 
of the nonzero parameter V. For simplicity, we will take the values of parameters as 

6=1, mi = m 2 = 1, and k 2 — 1 (8) 

in subsequent discussions. 

If the two oscillators, with fixed total energy, are separated larger enough, apparently, 
there might be no collision at all. Hence, there should exist an energy threshold E c . 
For a given energy E, the system is integrable as E < E c , and the system is in general 
nonintegrable as E > E c . The energy threshold corresponds to the case of xi = x 2 = x and 
Pi — Pi — 0, at which the collision just takes place, without exchanging the momentum. 
We know that the energy in this case can be written as a function of colliding position x 
and the separation parameter 6=1: 

E(x, 6 = 1) = hi(x + l) 2 + h 2 (x - l) 2 . (9) 

For the parameter 6=1 fixed, E(x, 1) has a minimum value with respect to x — x c . We 
thus have 

ki(x c + 1) + k 2 (x c - 1) = 0. (10) 
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x c can be obtained 



(h - fci) 

h + k 2 



(11) 



Therefore, substituting x c into eq.(9), we have 



E c = E(x c , 1) 



2h 



< 2, 



(12) 



fci + 1 



where the conditions in eq.(8) are used. The system will be nonintegrable if the initial 
energy E > 2 for any k±. This result is in agreement with the conclusion given in Ref. [13]. 
We will only consider the case for E > 2. 

To investigate the classical behaviors of the system, it is needed to obtain the corre- 
sponding Poincare sections. To realize this, we have to integrate out eqs. (7) numerically. 
Here it should be noted that the ratio of frequencies of the two oscillators is related to the 
ratio ki/k 2 by 



One may infer that the ratio uj\/uj 2 is rational or irrational, will give rise to different dynamic 
behaviors, or equivalently, different shapes of the Poincare sections. For this purpose, let us 
now present our numerical results. By integrating out eqs. (7) numerically, we have obtained 
the Poincare sections for k 1 /k 2 = 3, 4, 8, 9 with E = 2 and 5, and k\jk 2 = 3, 4 with E — 10 
and 30, respectively. In the calculation, we have kept the energy of the system to at least 
six effective digits. The Poincare sections for k\jk 2 = 3,4 and E = 2 are shown in Fig. 
1(a) and (b). Notice that in the calculation we have taken 60 different initial values, in 
the same region, for both ki/k 2 = 3 and 4. It is clear that the former shows some chaotic 
behaviors, while the latter shows the quasiperiodic behaviors, at least for the present initial 
values. This seems to exhibit that the ratio of ui/u 2 is integer or irrational, indeed results 
in different dynamic behaviors. However, further investigation shows that when the energy 
becomes larger, the differences between the integer ratios and the irrational ratios are vague. 
To show this, the Poincare sections for ki/k 2 = 3,4 and E = 5 are given in Fig.l(c) and 
(d) respectively. Unlike the cases of E — 2, the differences between the Poincare sections 




(13) 



5 



are ambiguous for E = 5, i.e., the ratio of uj\/uj 2 is integer or irrational, may not affect 
significantly the dynamic behaviors of the two colliding oscillators as the energy becomes 
larger. To verify above statement, we have still calculated the Poincare sections for the other 
values of k\jk 2 and E. For instance, the Poincare sections for k±/k 2 = 8,9 and E = 2, 5 
are obtained, as shown in Fig.2. One may observe that it shows the similar behaviors as 
k\jki = 3,4. The Poincare sections for ki/k 2 = 3,4 and E = 10,30 are also calculated, 
but they show nothing new except similar behaviors as E = 5. All above analyses mean 
that the ratio of frequencies of the two oscillators is integer or irrational, will influence the 
dynamic behaviors of the system drastically only when the energy is near the threshold E c 
(notice that E c = 1.6 for hi — 4 and 1.8 for k x = 9), i.e., the integer ratio exhibits some 
regular orbits, and the irrational ratio presents some chaotic orbits, as the energy is near 
E c . Whereas the integer or irrational ratio Ui/u 2 does not play essential role in the classical 
dynamic behaviors of the system, as the energy is far away from the threshold E c , i.e., the 
system can exhibit some chaotic behaviors either for the irrational ratio or for the integer 
ratio in the far as we have studied. It should be pointed out that in above discussion 

we usually start from the same initial conditions for given energy E. 

In addition to, to study the chaotic behaviors of the system in detail, we have also cal- 
culated the largest Lyapunov character exponent (LCE) [14] for different k±/k 2 and different 
energies away from the threshold E c (E > 2) in some chaotic regions. The result, as shown 
in Fig. 3, shows that the largest LCE depends on k 1 /k 2 linearly 

A max = P^, (14) 
k 2 

with P a constant, and is independent of E. This suggests that the classical chaotic behaviors 
of the system are mainly controlled by the ratio of frequencies of the two oscillators, as the 
initial energy is away from the threshold. 
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III. THE QUANTUM CASE 

Now let us study the quantum mechanical behaviors of the two colliding harmonic oscil- 
lators. The Hamiltonian of the system is 

H = H + H', (15) 

where H = H± + H 2 , Xi, Pi {i = 1, 2) satisfy the canonical commuting relations, and eq.(7) 
should be replaced by its corresponding quantum partner. Evidently, H can be exactly 
solved: 

H \nm > = E nrn \nm >, (16) 



with the eigenfunction 



where 



\nm > = (p n (xi + l)(f m (x 2 - 1), (17) 



<pi(xj) = N t exp( — ^-x^H^ajXj), 



and Hi(x) is the Hermite polynomials. {\nm >} thus constitutes a complete ortho normal 
set. The eigenvalue E nm is 

E nm — Tiw\(n + -) + hio-iim + -), n,m = integer (18) 

where lo\ = \fk[ and u 2 = \fk 2 - However, we must solve, at present, the following 
Schrodinger's equation 

Hip(x 1 ,x 2 ) = Eif}(xi,x 2 ), (19) 

with E the eigen energy of the system. To achieve this, we can write i[)(xi,x 2 ) as the linear 
combination of \nm > : 

N 

ip(x 1 ,x 2 ) =J2 a nm\nm>. (20) 
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Because of the orthogonal property of \nm >, we have 

N 

{H nm ,ik — E5 n i5 mk )a nm = 0, 



(21) 



n,m 



where 



H nm ,ik =< nm\H'\lk > +E nm 5 n i5 m k, 



(22) 



and 



5; 



0, if n 7^ m 



nm 



1, if n = m 



Since it is difficult to obtain the analytic solution of eq.(21), we solve it numerically. The 
integration in eq.(22) can be carried out by means of the properties of ^-function and the 
Gauss-Hermite integration methods. Note that the matrix which consists of H nm ^ is real 
symmetric. To find the energy spectra, the eigen equation must be solved. The maximum 
number of energy levels we have solved is 961. We find that the statistical properties of 
961 levels are qualitatively same as those of 200 levels. In the following, we will give the 
numerical results in detail, in order to show our consequences clearly. 



The distribution P(S) of the energy level spacings is very interesting. We take ft 2 = 0.2. 
For ki/k 2 = 3 (The ratio of frequencies is thus an irrational number), the distribution of 
the energy level spacings obeys the Wigner distribution for different interaction parameter 
V, and exhibits the energy repulsion[15]. However, for k±/k 2 = 4 (The ratio of frequen- 
cies is an integer), the distribution of the energy level spacings shows the transition from 
Poisson distribution to Wigner distribution if V varies from 0.0001 to 0.1. Fig.4 gives a 
comparison for the above two cases. Notice that all figures in Fig.4 are given for 200 levels. 
For k 1 /k 2 = 8,9, the similar phenomena occur. These results can be compared with the 
case of two dimensional harmonic oscillator where the P(S) does not exist if the frequen- 
cies are commensurable, and is peaked about a nonzero value of S if the frequencies are 
incommensurable [16]. 



A. P(S) versus S 
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B. \tp(x 1 ,x 2 )\ 2 

The density of probability of the system at (xi,x 2 ) plane is given by \i/j(x 1 , x 2 )\ 2 . Here 
Ti 2 is fixed to be 0.2. The surface plots and the related contour plots are given in Fig.5. 
The system shows some localization. It seems that the density of probability is more local 
for the case of k±/k 2 — 4 than for the case of ki/k 2 = 3. The scars in Fig. 5 almost lie in 
the region which is very near the accessible region for the classical orbits, but there is no 
clear relation between the density of the probability and the classical orbits, like the case of 
billiard systems [17]. For different V and k±/k 2: both the eigenvalues and the eigenvectors 
are so different that it is difficulty to compare the density of the probabilities for different 
parameters. The density of the probabilities has no explicit transition when the energy (i.e., 
the eigenvalues) crosses the classical chaotic energy threshold. 

C. E(V) versus V 

The relation between the eigenvalues and the interaction constant V is also studied, h 2 
is still fixed to be 0.2. We find that the eigenvalues are only sensitive to very small values of 
V. When V is larger than 1, the eigenvalues do not change with respect to the interaction 
parameter V. When V approaches zero, there is no degenerate energy level for k±/k 2 = 3, 
while there exist clearly the degenerate energy levels for k\jk 2 — 4, as shown in Fig. 6. In Fig. 
6, except V = 0, there is no energy level crossing, even though several lines are very close. 
For ki/k 2 = 8,9, the similar phenomena occur. We have also calculated the corresponding 
energy curvature [16, 17] (— ^|) numerically where the average is for 50 different levels, as 
shown in Fig. 6. The curvature seems to be independent of k\jk 2 . This result can be fitted 
by (1 ?v2)2 for small V, as given in [18]. 
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D. E(ki/k 2 ) versus k\jk 2 

Now let us investigate the relation between the eigenvalues and the ratio k\jk 2 . h 2 is yet 
fixed to be 0.2. V is fixed to be 1. The energy repulsion is very clear, as shown in Fig. 7, 
where k\ takes the values from 1 to 25. In Fig. 7, there is no energy crossing in the given 
region as well. We have calculated the curvature (— ) (k 2 = 1) numerically where the 
average is taken for 50 different levels. The result can also be fitted by jj^zjz- 

E. E(h) versus fi 

The relation between the eigenvalues and the quantization constant H is studied also. In 
this case, the parameter V takes two values: 0.001 and 1.0. The parameter k\jk 2 also takes 
two values: 3 and 4. For ki/k 2 = 3, the energy as the function of Ti is almost independent 
of the values of the parameter V. However, for k\jk 2 = 4, the energy as the function of h 
is quite different for different values of the parameter V, as shown in Fig.8. For k\jk 2 = 4, 
when V is small, the degenerate energy level occurs as H is large, and when V is large, the 
degenerate energy level occurs as H is small. We still can not fully understand this surprising 
phenomena for the moment. As expected, the semiclassical quantization constant H plays a 
crucial role in discussions of the transition from quantum systems to classical systems. 



IV. DISCUSSION AND CONCLUSION 

We have systematically studied both the classical and quantum chaotic behaviors of two 
colliding harmonic oscillators. In the classical case, there exists an energy threshold above 
which the system can be chaotic. The system is conserved, and the classical behaviors 
of the system belongs to the KAM class. Moreover, the ratio of the frequencies of two 
harmonic oscillators affects the Poincare sections significantly, as the initial energy is near 
the threshold. The largest LCE depends linearly on the ratio ki/k 2 of the two oscillators, and 
is independent of the initial energy in some chaotic regions, as the energy is away from the 
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threshold. This implies that the chaotic behaviors of the system are mainly controlled by the 
ratio ki/k 2 , not by E in the situation. In the quantum case, the ratio of frequencies of the two 
harmonic oscillators is also important. Like in the two dimensional harmonic oscillators, the 
incommensurable frequencies play a special role. When we deal with the distribution of the 
energy level spacings, the integer ratio of frequencies shows the transition from the Poisson 
distribution to the Wigner distribution with increasing V. If the interaction parameter 
V is not very small, the properties of the system are independent of V. The relation 
between the energy levels and the ratio of frequencies of the two harmonic oscillators, and 
the semi-classical quantization constant H is investigated, respectively. The curvature of the 
energy levels with respect to V or k±/k 2 shows universal properties. As one notes that, the 
parameter V plays no role in the classical case, because the potential barrier at x\ = x 2 is 
infinite, in which makes V playing no role. While in quantum case, there are possibilities 
for a particle transmision over an 5-function-type potential barrier, and V should survive. 
Our present results just confirm this fact. Because of the numerical difficulty, the total 
behaviors of the system in the processes of H approaching zero have not been obtained. 
Nevertheless, the partial results have already shown some very interesting phenomena. One 
may observe that such an apparently simple system indeed exhibits rather complicated 
dynamical behaviors, and a few phenomena can not be sensibly explained for the moment. 
This system may yet contain rich exotic behaviors, and deserves to further study. How to 
understand the transition from a microscopic system, through a mesoscopic system, to a 
macroscopic system, is still a fascinating topic. 

In summary, for some values of the initial energy near E c , we have found that the 
ratio of frequencies of the two oscillators affects the Poincare sections significantly. The 
largest Lyapunov character exponent depends linearly on the ratio of frequencies of the two 
oscillators away from the energy threshold in some chaotic regions, which shows that the 
chaotic behaviors of the system are mainly determined by the ratio. In the quantum case, 
for certain parameters, the distribution of the energy level spacings also varies with the 
ratio of frequencies of the two oscillators. The relation between the energy spectra and the 
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ratio of frequencies of the two oscillators, the interaction constant, and the semi-classical 
quantization constant, is also investigated respectively, and some interesting phenomena are 
presented. 
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Figure Captions 

Fig.l. The Poincare sections for k-i — 3,4 
Fig. 2. The Poincare sections for k\ = 8, 9 

Fig. 3. The largest Lyapunov character exponent versus ki/k 2 (k 2 = 1) for different 
gies E > 2 

Fig. 4. The distribution of level spacings for different V and k\. 

Fig. 5. \ip(xi 1 X2)\ 2 versus x±, x 2 (The wave functions are not normalized.) 

Fig. 6. The Energy E and the curvature (— ^) versus the interaction constant V 

Fig. 7. The energy E and the curvature (— ) versus k\ 

Fig. 8. The energy E(h) versus h 
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